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Abstract 

Motivated by the inverse problem for the Lemaitre-Tolman-Bondi dust solution, in 
which problem the luminosity distance function Dl{z) is taken as an input to select 
a specific model, we compute the function Dl{z) of the LTB solution up to the third 
^SJ . order of z. To perform the otherwise cumbersome computation, we introduce a new 

' convenient form of the LTB solution, in which the solution is explicit and unified. 

With this form of the LTB solution we obtain the luminosity distance function with 
' the full generality. We in particular find that the function exactly coincides with that 

. of a homogeneous and isotropic dust solution up to second order, if we demand that 

' the solution be regular at the center. 

■ 

o 

^ ' Modern observations show that the Universe is presently in an accelerating 

phase and dominated by dark energy with negative pressure [1-5]. It is re- 
markable that there seems no obvious contradiction among the observations, 
but the true nature of dark energy is a great mystery. On the other hand, it 
' is also true that this recognition is a result of the strict use of a Friedmann- 

Lemaitre-Robertson- Walker (FLRW) homogeneous and isotropic model. If in- 
homogeneities are properly taken into account it might be possible to explain 
the observations without introducing dark energy. Recently, this possibility has 
renewed interest in inhomogeneous cosniological models, especially in the so- 
called Lemaitre-Tolman-Bondi (LTB) solution [6-8], which represents a spheri- 
cally symmetric dust-filled universe. Because of its simplicity this solution has 
been considered to be most useful to evaluate the effect of inhomogeneities in 
the observables like the luminosity distance-redshift relation. 

Although this solution is simple, it has great flexibility. Various models 
have been proposed that are consistent with the observations [9-15]. But, what 
is the best configuration that can explain the observations of, e.g., luminosity 
distance redshift relation for type la supernovae and still be consistent with other 
observations? To respond to this question one needs a systematic approach. 

*E-mail: tanimotoOgravity .phys .nagoya-u. ac . jp 
^ E-mail: nambuOgravity . phys . nagoy a-u .ac.jp 
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Our primary focus is the "inverse problem" approach, in which one takes a 
luminosity distance function Dl{z) as an input to select a specific LTB model. 
Mustapha et.al. [16] argued that the free functions in the LTB solution can 
be chosen to match any given luminosity distance function Dl{z) and source 
evolution function m{z). Celerier [9] performed an expansion of Dl{z) for the 
LTB solution of parabolic type to fit the ACDM model with = 0.7 and 
~ 0.3, and argued that the model can explain the SNIa observation at least 
for z < 1. Vanderveld et.al. [17] however showed that such a fitting with an 
accelerating model is only possible at the cost of occurrence of weak singularity 
at the center. 

The main purpose of this paper is to find a Taylor expansion of Dl{z) for 
the LTB solution in the most general form, not restricted to the parabolic type, 
with considerations for the condition to avoid the weak singularity. Since the 
differential equation for the inverse problem is in general singular at the origin 
z = it is imperative to prepare a solution there using a different method. If 
we have the expansion of Dl{z) for the general LTB solution, we can easily find 
this solution by comparing the given Dl{z) and that of the LTB solution order 
by order. 

To perform the computation we present a new way of representing the LTB 
solution. Although this is not our main purpose we would like to stress that this 
new representation would be very useful in various computations concerning the 
LTB solution. This solution is usually expressed with parametric functions, as 
is the FLRW dust solution. This parametric character of the solution often 
makes our analysis and considerations complicated and non-transparent. More- 
over those functions change their forms depending on whether the solution is of 
spherical type, parabolic type, or hyperbolic type. This split-into-cases charac- 
ter is also a drawback of the conventional expression of the solution. The new 
expression of the solution dissolves all these unwanted characters. It will play 
an essential role to achieve the complicated computations needed to compute 
Dl{z) in the fully general setting. 

The structure of this paper is as follows. In section [2] we introduce the new 
unified form of the LTB solution. In section [3] we perform the expansion of 
Dl{z) using the unified expression. In section |3] wc study how different the 
function Dl{z) for the LTB solution is from that of the FLRW dust solution. 
Section [5] is devoted to conclusion. Appendix 1X1 presents the result without the 
regularity condition at the center. 



2 The LTB solution in a unified new form 

Let us first gather the conventional forms of the solution. The LTB metric is 
written in the form 

1 -I- 2E{r) y ^ ' y ^ w 

where primes denote derivatives with respect to the radial coordinate r, and 
E{r) is a free function called the "energy function" . The conventional way of 
expressing the areal radius function R(t,r) depends on the sign of E{r) and is 
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parametric: For E{r) > 0, we have 



R{t,r) (2) 
M(r) 

t^tBir) ^ ^^^^^^^3^^ (sinh??-r;), (3) 

where M{r) and tsir) are free functions caUed, respectively, the "mass function" 
and the "big bang function" . For E{r) < 0, we have 

R{t,r) = ^^il~co,rj), (4) 
M(r) ^ 

t~t,ir) = ^-^^-l^i,^s.n,). (5) 
Finally, for E{r) — 0, we have 

R{t,r) = (^^y Miry/^t^tB{r)f^'. (6) 

(In this case the function R(t, r) is explicit in terms of t and r.) 

Now, observe that in the E > case, Eq.Q shows that 77 can be regarded 
as a function of 2E{{t — tB)/M)'^/^. Therefore from Eq.([2]) the function R can 
be written in the form 



using a certain function X{x). To make the E ^ Q limit apparently regular, we 
then factor out x from the function X{x) and make X{x) = Q^/^xS{—Q^'^/^x)] 



R{t,r) = m{r)f''{t-tB{r)f''S (-2E{r) (^^^) 



2/3N 



(8) 



(The numerical factors inserted are our convention.) 

The same observation is also applicable to both E <Q and E = Q cases, and 
as a result, we find that the form ^ may provide a desirable form of R{t,r), 
which is explicit in terms of t and r and requires no separate considerations 
depending on the sign of E{r). There is, however, still a remaining task, which 
is to confirm that the function S{x) is smooth at a; = 0. Otherwise, this form 
would be superficial. 

To this, note that the function S{x) can be expressed in parametric forms. 
For a; < (corresponding to i? > 0), we have 

' cosh^/^-l ^_ -1 r-: r-^x2/3 



'^(-)=6V^ (sinhV^-W/3 ' --^(^-W-C-V-C)^^^. (9) 
For a; > (i? < 0), we have 

*>^6^^7fs^> .^5^(v^-.inV?)^'^ (10) 
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For a; = 0, we have 



m = ( 4 



1/3 



(11) 



The parameter C takes positive and negative values in accordance with the sign 
of X. (The expHcit correspondence to rj is given by C = 7]^ for E < 0, and 
C = — for E > 0.) We then find that S and x are expanded in powers of C in 
the same form for all signs of x, since by direct computation we can immediately 
confirm the following common expression: 



Six) = 



1 



1 - LC + 

1 4!*' ^ 



61/3 



L 

5! 



2/3 ' 



62/3 



5! 



2/3 



(12) 



Proceeding the series expansions, 



S{x) = 



1/3 



1 - 



20 



1680 

3 



o(C') 



30^ 25200 +^^^ 



(13) 



Since both S and x are smooth (in fact, analytic) functions of C, in particular 
in the neighborhood of C = 0, this shows that S(x) is smooth at cc = 0, as 
desired. (The smoothness ol S{x) at x 7^ is apparent.) Figure [T] shows a plot 
of the function S(x). S{x) is a non-negative, monotonic function, defined for 
x<xc = (V3)'/^- 



2.5r 




Figure 1: Plots of the function S{x) (solid line), the function \>q(x) — \J~x 
(dashed line), and the function Sc(x) = \Jxc — x (gray line). Although the 
parametric expression of S{x) takes different forms for a; > 0, a: = 0, and a; < 0, 
this function is smooth at the dividing point a; = 0. 

The expression (jH]) therefore does provide a desired form of the LTB solution, 
together with the conventional metric form ([T]). This expression will turn out 
to be extremely useful in performing most of the computations concerning the 
LTB solution, including the expansion of Dl{z) in the next section. In the rest 
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of this section we briefly discuss some useful properties concerning the function 

First, remember that R(t, r) satisfies the foflowing generahzed 'Friedmann' 
equation: 

R^^ — + 2E, (14) 

where dots denote derivatives with respect to the proper time t. Substituting 
Eq.Q we immediately have the first order ordinary differential equation (ODE) 
for S{x): 

\{S{x)+xS'{x)f +?>x-^^Q, (15) 

where S' = dS{x)/dx. (We understand that primes attached to S always stand 
for derivatives with respect to its single argument, not to the radial coordinate 
r). 

This ODE has a peculiar feature observed at a: = 0; at this point the term 
xS'{x) vanishes and therefore the equation degenerates into an algebraic equa- 
tion that constrains 5(0), provided that |5'(0)| < oo. As a result we obtain 
5(0) — (3/4)^/^ in accordance with Eq. fTTj) . The finiteness condition of S'{x) 
at a; = implies that the line y = S{x) in the x-y plane should intersect the 
X = axis transversely. The above feature therefore tells us that the function 
S{x) is the unique 'transversal' solution of the ODE il5\) . It may be useful to 
use this characterization to define S{x), instead of using the explicit parametric 
expressions ©, (dH]), and (fTTjl . 

It may be worth commenting that we should not think of the ODE (fT5|) 
as a usual 'evolution equation,' because of the following two reasons. First, 
as discussed previously, this ODE has no freedom to choose an initial value. 
Note that the solution of the generalized Friedmann equation does possess 
freedom to choose an initial data, which is the function tsir). This function is 
however already incorporated in the form ([H]). The resulting equation therefore 
cannnot contain further freedom to choose a solution. Second, the gradient of 
the variable x{t, r) for the ODE is in general not timelike. In other words, the 
contours of x{t,r) are not always spacelike. This is most noticeable when E(r) 
crosses zero at some spatial points. Figure [5] shows an example where E(r) 
crosses zero at one point (r = 1). 

Note that the form ([8]) already contains the parabolic solution as one factor, 
with the remaining factor being S{x). We can think of \x\ as a parameter that 
measures the deviation from the parabolic evolution. It is obvious that x — 
corresponds to both big bang time t = t_B(r) and comoving points E{r) = 0. 
In the above example, these two regions are shown as the bold lines forming an 
inverted 'T'-shape (Figl2]). Let us consider a neighborhood of this 'T'-shaped 
region such that for a small positive e, |a;| < e. Then we can say that the time 
evolution of the points contained in this spacetime region is close to that of 
the parabolic solution. Actually, for small t — tsir) this is a well-known fact. 
(Both elliptic and hyperbolic solutions asymptotically approach the parabolic 
type evolution {t — tsir))^^^ near the big bang time.) It is also apparent that 
if |i?(r)| is small enough the evolution of those points is close to that of the 
parabolic solution (for a finite time- interval) . The variable x is therefore a 
deviation parameter, rather than a evolution parameter. 
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Figure 2: Typical contours of x{t, r) when E{r) crosses zero. The decreasing 
bold curve corresponds to the Big Bang singularity {t — tsir)), while the vertical 
bold line to E{r) = (this is not a singularity). Both bold lines correspond to 
X = 0. Contours are drawn for t > tsir). [This example is intended to be helpful 
for a general discussion of the properties of S{x), and not to imply a setting of the 
next section. The specific choice for this example is E{r) — —r'^ {I /2—do,4{r—l)) , 
M{r) — r^, and = e^''~ . The function 9c{r) is a step function with the 

transition interval [— cr, a] with values 0(j{r) = l(r > a), and 0(r < — u).] 

Now, let us return to the discussion of useful properties of S(x). Remem- 
ber that Eq. (fT4|) is an integral of the following second order partial differential 
equation: 

2RR + if -2E = 0. (16) 
From this we have the second order ODE for S{x): 

x{2S{x)S"ix) + S'\x)) + 5S{x)S'{x) + ^ = 0. (17) 

For X ^ 0, this equation is useful in case we want to eliminate higher derivatives 
of S{x) than first order. 

When X = 0, Eq. (fT7|) itself does not determine S"{Q), due to the multiplied 
factor X. However, we can regard this equation as an equation that determines 
S'{0) from 5(0), and find S'{0) = -9/(205(0)). Similarly, taking derivatives 
of the equation we can successively determine the values of arbitrarily higher 
derivatives of S{x) at a; = 0, in terms of 5(0). One of the useful applications of 
this property is to obtain a series expansion of S{x) about a; = 0. Since 5(0) is 
given in Eg. pip (or determined from Eq. (|15p ). we have 

As discussed previously, this expansion can provide a 'parabolic approximation' 
of the solution. 
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We may be interested in the significance of the general solutions of the ODE 
PT)) . To avoid confusions let us top tilde, like S{x), to denote a solution of the 
ODE which does not necessarily coincide with S{x). First, the above observation 
immediately tells us that a solution of the ODE which intersects the x = axis 
transversely is specified only by one parameter s = 5(0), since S'{x) is not 
free at a; = 0. We can easily confirm that all those transversal solutions are 
generated from S{x) by rescaling: 

S^ix) = aS (^) , (19) 

where a = s/S{0) = (4/3)^/'^s. This rescaled function is also a solution of the 
first order ODE (fT5|l with the modification 1/5 ^ l/{a^S). These rescaling of 
the function and modification of the equation just correspond to the rescaling 
M(r) — > a^M{r), and therefore although the rescaled function Sa{x) does gen- 
erate a dust solution, it is equivalent to the one by the original S{x). Thus, we 
may not be interested in the general transversal solution Sa{x). 

Non-transversal solutions of the ODE (fT7|) in general have two-parameters. 
An interesting one-parameter special solution is 

V;3(x) = V^-^, (20) 
x 

where /3 is an arbitrary constant parameter. Vf3{x) approaches S{x) asx ^ —oo: 

Vp{x)^S{x) (x^^oo). (21) 

In particular, Vo(a;) = \/—x approaches S{x) from below (FiglT|), and is often 
useful for various estimates of S{x) for large —x. We remark that the function 
Vf3{x) satisfies the first order ODE (fT5|l if the 1/S{x) term in the equation is 
neglected. This in effect corresponds to taking the limit M(r) 0, and therefore 
the metric generated by this function (through Eq.([8]) with S{x) replaced by 
Vfs^x)) represents a vacuum solution. In fact, by a direct computation of the 
Riemann tensor we find that it all represents the Minkowski solution. 

Finally, it is worth mentioning that, slightly modifying Vo(a;), the function 



Scix) = y/xc ~ X, (22) 

which is no longer a solution of any of the ODEs, provides a good approximatiorQ 
oiS{x) for all domain of a; < xc- The function Scix) asymptotically approaches 
S{x) from above as x — oo (Fig[T|). An elementary estimate actually confirms 
the inequality 

Six) < Scix), (23) 

where the equality holds for x — xc- Also, another elementary estimate estab- 
lishes 

xS'ix) < ^ (24) 

for the first derivative. 

^It is interesting to note that in the special case of the FLRW dust solution, the approx- 
imated solution with Sc{x) corresponds to the "renormalized solution" [18] of a solution 
obtained in the long wavelength approximation of Einstein's equation. 
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3 Expansion of Dl{z) 



Motivated by the inverse problem, we in this section perform the expansion of 
the function Dl{z), the luminosity distance as a function of the redshift, for the 
LTB solution. 

The luminosity distance for the LTB solution is simply given by [19,20] 

Dl = {1 + z)^R, (25) 

where the areal radius function R = R{t{z),r{z)) must be evaluated along the 
light ray emitted from the light source and caught by the central observer at 
{t, r) — (to, 0). With this formula we can expand Dl{z) if we expand R. 

Before proceeding the expansion, we note that the LTB solution has the 
gauge freedom of choosing the radial coordinate r. In fact, it is easy to see that 
coordinate transformation r ^ f = f{r) retains the characteristic form of the 
solution if we redefine the free functions E{r),tB{r), and M{r) suitably. The 
significance of this freedom is that it enables us to fix the function M(r) for 
many cases. 

One of the popular gauge choices is (e.g., [9, 14, 17]) to take M{r) = Mor^ 
with Mq being a positive constant. We call this choice the FLRW gauge, since 
the standard form of the FLRW dust solution is realized with this choice with 
the other functions being E{r) = — (fc/2)r^, and tsir) — 0. {k is the curvature 
constant.) 

Another excellent choice is [16,21] to choose the radial coordinate r so that 
the light rays coming into the observer are simply expressed by 

t = ta-r. (26) 

This is equivalent to the condition that dr/dt = —1 must hold along those light 
rays, and therefore from the metric ([1]) it is equivalent to imposing 



R'{to-r,r) 



= 1, (27) 



where i?'(to — = {dR{t,r) / dr)\t=to-r- This gauge choice determines the 
function M{r) only simultaneously with E{r) and tsir). We call this choice of 
r or the resulting choice of M{r) the light cone gauge. 

We in the following adopt the light cone gauge, which, together with the 
formula ([H]), makes our procedure of expansion systematic and straightforward. 
In fact, with the condition (|26p we can firstly expand R in powers of r in a 
straightforward manner, putting R ~ R{to — r, r) in the formula ([5]). To convert 
the result into expansion in terms of z, we need to find r(z) in powers of z. This 
is possible from the formula [9] 



^ = + ^ . (28) 

dz {1 + z)R'{tQ~r,r) 

This would complete our expansion, but we still need to consider an extra 
condition, which is the regularity at r = 0. As pointed out in [17,22], if the 
first derivative with respect to r of the matter density p(t, r) does not vanish 
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at r = 0, the spacetime has a weak singularity there. We wish to ehminate this 
undesirable feature; therefore impose 



p'(i,0)-0. (29) 

The expansion with this condition imposed will give the final form of our result. 

We are now in a position to proceed the expansion, following the procedure 
outlined above. We wish to expand Dl{z) up to the third order, which leads us 
to expand the free functions as follows: 

tB{r)=b,r + ^r^+Oir^), (30) 

where the differential coefficients a = E'^'\0), b, = (0), and m, = M(*)(0) 
are constants, with 7713 > 0. Our task is then to express the expansion of Dl{z) 
in terms of these coefficients. 

It turns out to be convenient to, prior to the general computations, finish 
the computation for the first order of Dl[z) to determine the Hubble constant 

in comparison with the definition Dl(z) = z/Hq + O(z^). This is not very 
much elaborating, and we find 

■^o==:7r 1 + 2^0-^7 — ^ ' xo = -e2{ — . 31) 
3to V S[^o) ) V"^3/ 

(As remarked, prime attached to S(x) stands for the derivative with respect to 
cc, not with respect to r; S' (xq) = {dS[x) j dx)\x=xo^ 

The gauge condition ((27|) is equivalent to the vanishing of the function 



G(r) = v/1 + 2^(r) - Blit^ -r,r). (32) 

We expand this function up to second order and demand that the coefficients 
be equal to zero. First of all, from the 0th coefficient we have 

S{xo) =tg^^^m^^^^, (33) 
and together with Ea. ((5T|l . we also have 

S'M = ^^"^^tH'' . (34) 

In the following we will use Eqs. ([55|) and ([M]) to eliminate S{xo) and S'{xo) from 
various equations. We remark that as we will see, Eq. ((M| is valid not only for 
the case 62 7^ but also for the 62 = case, if we understand that an appropriate 
limit is taken. Note also that Eqs. ([55|) and are not independent, since so 
are not S{x) and S'{x). Substituting the equations into Eq. (|15p we obtain the 
relation 

772 = 1 - T772- (35) 
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Since > 0, this also means 

< 1. (36) 

To take the 62 — > hmit in various equations that wiU appear in the rest of 
our computation, it is useful to have an expression of Hota in powers of 62 /Hq. 
First, from Eqs.JST]) and HH) we have 

^ 3 /3\^^^ 117 



On the other hand, from Eq. ((55l) (and the definition of xq (pij) ') we have 

^0 = -T72 ^ -TT72 ■ (38) 



ff2 Vl-e2/7?o' 

Substituting this into Eq.(P7)). we obtain an equation for iJo^o- We can solve 
this and find 



This is the expression we wanted. In particular, this implies as 62 — ^ 0, 

2 — 3-ffo^o 2 

2-3i?oio 2\ 1 6 



e2/i/2 5; e2/i/o' 35' 



(40) 
(41) 



With the first of these we can confirm that Eq. provides the right value for 
the limit 62 0, and therefore Eq. ([M)l is, as mentioned, valid also for 62 — 0. 
Returning to the gauge condition ((5^ , from its first order coefficient we have 



1/3 

4m3 4627713' S'{xo) 
7774 = — ttI = GHqUIj. (42) 

Similarly, from the second order we have 

= ? (l5Hj + 5| + 5^^-%^ {UHS + 12H„., + 2e. - Si)) . 

(The second derivative S"{xo) appears in the coefficient, which we eliminate 
using Eq. pT]) .) In particular, for 62 = we have 

7775 = 7773 (^lOHoiAHo + 362) - . (44) 



In the following we will eliminate 7774 and 7715 with Eqs. (|42p and (|43p . 
To find the regularity condition at r = we expand 

M'(r) 

4np(t, r) ^ — ^ (45) 

' ' R^{t,r)R'{t,r) ^ ' 
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in powers of r up to first order with t fixed: 
1 



47rp(i, r) 



S^{x{t)) 
( 46i 



1 
2fi 



46ie2 S'{x{t)) ^ (2637713 - e2TO4) S'{x{t)) 
iml'h^/^ S{x{t)) 



,3t3 2,ml'hV^ S{x{t)) 



r (46) 



+ 0(r2) 



where 



The condition 



then implies 

bi = 0, 63 



eat 



2/3 



2/3 



627774 
27ri3 



3iJoe2- 



(47) 



(48) 



In the last equality we have used Eq. ((l^ . (If we had employed the FLRW 
gauge, which implies 7714 = 0, we would have had 63 = and 61 — 0, which is 
consistent with the condition derived in [17].) 

From Eq. (|^5)) we can immediately find the energy density at the central 
observer. Using Eg. ([55)1 we have 



87rp(to, 0) = 7773 
Therefore it may be natural to write 



3Hl 



-o 

^0 



(49) 



(50) 



where 57j\/.o and flk,o are density parameters for the central observer. Then we 
can interpret Eg. ([55)1 as the usual relation for the dust FLRW model: 



n 



M.O 



(51) 



The final task before presenting our main result is to obtain r(z) up to third 
order. Note that Eq. (P5| is integrable and we find 



I + z = exp 



y/TTmy) 



dr'. 



(52) 



From this equation, we have 



1 



Ho iHo 
1 

1 — HqIq 



5- 



62 
-"0 



48i/o 



109 



35e| + 364 46^ + 64 



^0 



13 



1262 116^+64 126262 26^-64 



Ho 



Hle2 



H^ 



H^ 



z3 + 0(z4). 
(53) 



In particular, for 62 = we have 



Kz) 



1 /35 96, 



Ho AHo 



6H0 



64 



4 2Ho lOH^ 



z^ + 0{z^). (54) 
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The regularity conditions have been imposed in these equations, as well as 
the gauge conditions (|42p and 

We are now ready to compute the expansion of R in the light cone gauge in 
powers of r with the regularity at r = 0, and convert it to in powers of z using 
the formula ([53|) or (I54p . As mentioned, the expansion of R — R{to — r,r) is 
straightforward if we write R using the function S{x) as in Eq.®. Moreover, the 
use of Eas. ([55|) and ([M)) . together with the ODE ([T7)) . enables us to eliminate 
the function S{x) and its derivatives to simplify the equations. 

The result of the expansion of R is in such a simple form that does not 
depend on the higher differential coefficients 64 and 62: 



R{to ~r,r)=r - —r- - ^ 



1 - 



62 



0{7 



(55) 



Substituting Eq. (l551) and taking into account Eq. ipS)) . we have the final form of 
regular Dl{z), which is 



1 + ^ 

e2/gg ^3^262 



12&, 



1 



1 — Hoto 



Hoto 



Ho 
2 ~ SHptp 
e2/Hl 



/ 4e2 



64 

H'^ 



(56) 



(This looks singular when 1 — Hpto = 0, but the use of Eqs.((3T|l and (|24|) 
immediately shows 1 — Hoto > 0.) In particular, for 62 = 0: 



Dl{z) = — + 



1 



Ho 4Ho 



8H0 



662 
Ho 



264 

15H^ 



z3 + 0(z4). 



(57) 



It may be of some interest to have Dl{z) without the regularity conditions, 
which is presented in Appendix \K\ 

With the above result we can easily determine an LTB model so that its 
Dl{z) coincides with a given luminosity distance function £)^'"^(z) at least up 
to third order. The input function d'"^^\z) may be expanded as 



Dl \z) = hz + —z +-^^ + 
A comparison with Eq. (|56p immediately gives 



(58) 



Hq — —, 

62 _J2 



1, 



(59) 
(60) 



and 



1 



Hot 



Hoto 



2-3i^o<o^ 64 , 1 - 62/gg 1262 
1 — Hoto Ho 



e2/Hl 

462 
/2 



H^ 



1 - 62/gg 
1 — Hoto 



13 



262 
-"0 



(61) 



12 



Equation (j59p determines Hq, which constrains the parameters 62, m^, and to 
through Eq. (j3ip . Equation (|60p determines 62, which in turn determines ma 
through Eg. ((55)) . Then, the present time to for the central observer is imphcitly 
determined from the condition p3p with 62 and 7713. This way we can determine 
the parameters 62, ms, and Iq for given /i and I2. We can see that the right 
hand side of Ea. (j6ip is now a function of Ii, I2 and I^. This equation represents 
a constraint for 64 and 62, only one of which is determined from the other 
through this equation. This is a result of the fact that the LTB solution is 
highly degenerate in that it can represent inequivalent multiple models that 
give rise to the same Dl{z). If we determine one of 64 and 62 according to a 
separate consideration, we can determine the other. All the other parameters 
are determined from Eqs. (|^^ . and ([^5]) . We remark that this procedure is 
solvable only for /2//1 < 1 due to the inequality ([5^ . (See the next section for 
its significance.) 



4 Comparison with the FLRW Umit 



To understand the significance of the expansion ([Sd]) itself, it is profitable to 
consider how much different it is from the FLRW limit. As mentioned, one of 
the choices that realize the FLRW dust solution in the LTB solution is to take 

E{r)^-\r^^^r\ M{r) ^ M,r^ ^ ^^r\ tsir) ^ Q. (62) 

{tBif) can be any constant, for which we take zero.) However, this coordinate 
choice is different from the one we have been employing. A general gauge- 
independent condition to give the FLRW solution is obtained by eliminating 
the explicit r-dependences from the above expressions: 

(The constant factors have been determined so as to be valid for any M{r) — 
0{r^).) The functions that satisfy both these equations and gauge condition 
(P7|) provide the functions that realize the FLRW solution in the light cone 
gauge. Substituting the expansions ([30]) into the above equations we have 

63 = TT"^' ^4 = (-5777-4 + 2477737775) . (64) 

2r773 60777| ^ ' 

We must solve the four equations, the above two with Eqs. (|l^ and for the 
four variables 63, 64, 7r74, and m^. We have 



63 = 3i/oe2, 64 = (13i/o + 262)62, 
7774 = 6i/o?T73, 7775 = 5(8i/o + 62)7713, 



(65) 



which, together with hi — 0, give the FLRW limit for the differential coefficients. 
Among the above four equations, only the equation for £4 is essential; the equa- 
tion for 7774 is the same as the general gauge condition ()42p . the one for 7775 is 
equivalent to the general gauge condition (|43p with the equation for £4 imposed, 
and the one for 63 is the same as the regularity condition (|48p . (The regularity 
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at the center is therefore automatically guaranteed for the FLRW limit as it 
should be.) 

These limit conditions motivate us to put 

64 = {13H^ + 262)62 + H^e4, (66) 

where £4 is a dimensionless parameter which becomes zero in the FLRW limit. 
Substituting this into Ea. (j56p . we have 



48Fo(l - Hoto) 
+ 0(/), 



2-3Hoto \ ( 62 \ 1262 



(67) 



where the Dl{z) for the FLRW limit is 

A striking feature of our result is that it shows that the luminosity distance 
function Dl{z) Jot an LTB solution which is regular at the center exactly coin- 
cides with that of an FLRW dust solution, up to second order. 

In particular, if we, as in [23], define the deceleration parameter go in com- 
parison with the Dl{z) for a general FLRW model: 



Dl{z) = J^\^ + + 0{z') ) , (69) 



we have go > 0, as in the FLRW dust solution. In fact, comparing with Eg. ([56 
we have 

1 A 62 



2 V i?, 



90 = 1 - l4 > 0, (70) 



because of Eg . ([55]) . This reconfirms the claim in [17,24,25]. 



5 Conclusions 

We have first given a new way of expressing the LTB solution, which is explicit 
(i.e., not parametric) and reguires no separate considerations depending on the 
sign of the energy function E{r). This has been done using a "special" function 
S{x), which can be defined as the unigue 'transversal' solution of the first order 
ODE (|15p . Using this monotonic function, we can write the areal radius function 
R{t,r) for the LTB metric in the concise form ([8]). To simplify expressions 
involving higher derivatives of S{x) it is most useful to use the second order 

ODE (ini). 

Taking advantage of this concise expression, we have computed the luminos- 
ity distance function Dl{z) for the LTB solution up to third order of z. We 
have found that if wc impose the regularity condition at the center, the function 
degenerates into an FLRW dust case up to second order, and differences only 
appear from the third order. 
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The second order coincidence with the FLRW dust solution tells us that we 
cannot choose a set of LTB functions {E{r),tB{r), M{r)} so that the result- 
ing Dl{z) fits an FLRW model which shows an accelerating expansion like the 
ACDM model, as long as the regularity condition is imposed. This is however 
of course not to say that we cannot find an LTB model that explains the obser- 
vations, since for this it is not necessary to have an exact fitting of Dl{z) near 
the center with an accelerating FLRW model. 

Perhaps, the simplest way to guarantee the regularity and still have a good fit 
of Dl{z) with the observations is to choose the model so that it exactly coincides 
with an FLRW dust solution (of perhaps negative curvature) for z smaller than a 
certain small (but finite) value zi , and use the full flexibility of the LTB solution 
for z > zi to fit the observations. This approach however has the drawback that 
the Dl(z) chosen this way inevitably becomes a non-analytic function, since if it 
was analytic it would be that of the FLRW solution. To be sure, it is possible to 
approximate such a non-analytic function with an appropriate analytic function. 
For example, one might use tanha: to approximate a step function, but tanha; 
never becomes constant for large x, as opposed to the step function. Like this, 
an analytic Dl(z) chosen to approximate a Dl(z) that is endowed with the 
above property deviates from that of FLRW near the center. 

To summarize, if one wants Dl{z) to be analytic (or of any arbitrary form), 
the inverse problem at the center is nontrivial. Our result (|56p gives, as we have 
seen, a solution to this problem. A further study based on the present work is 
under progress. 
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A Dl{z) without regularity condition 

In this Appendix, we present the luminosity distance function for the LTB 
solution that is not necessarily regular at the center. To distinguish from the 
regular one, we write D^^^^\z) to denote this general function. Let Dl{z) be 

the regular part of d''£^^\z) as in Eq. ([56|) or ([67|) . and dL{z) be the difference 
from the regular part, i.e., 

Df\z)^DL{z) + dL{z). (71) 

The difference part dL{z) should vanish when the conditions (|48|) are satisfied. 
Motivated by the condition for 63 , we rewrite 63 using the dimensionless param- 
eter £3 defined by 

63 = 377062 + Hie3. (72) 
The parameter £3 vanishes when 63 satisfies the regularity condition. 
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Then, we have 

dLiz) = 

where 

A2 

and 



A, 



12Ho{l-Hotor ' 288Hoil ~ Hotor 
2 — 3Hoto 



Hoto + 



£3 + 6 ( 1- ) 61, 



(73) 
(74) 



A3=36(4-5i7oio) ( 1-;^ ) hi 



12 



2Hoto{S-7Hoto)-4. 



62 1(2 — Sffgio 



H^tliA-5Hoto)- 



2(2 — SHqIq 



{b- Hoh[2l-2QHoto)) 



hie 



2-3Hota , 2\ 10 - Hnto{22 - ISHoto) 
5 



£2 



62 



4(1 - Hoto)i 6 1-^9- llHoto - 2(3 - ^H^h)-^ h 



^0 



-"0 



e2/H^ 
A-2Hotoi3-AHoto) 



62 
^0 



£3 



(75) 



To take the hmit 62 — > of this expression, one may need both of Eqs. (|40|) and 
(|¥T|). As a resuh, for 62 = we have 



where 



and 



„ _ £3 36i 
^^=15 + —' 

7 £3 13 £3^ £361 5 61 9 61^ 



30 1575 2 4 4 
The deceleration parameter qo for the central observer is therefore 



1 



90 = 



1 - 



62 

^0 



1 



6(1 - Hoto) 



Hot 



2 — 3Hoto 
e2/H§ 



£3 + 6(1-^ 



In particular, for the 62 = case. 



(76) 

(77) 
(78) 

hi . 
(79) 

(80) 



which is apparently positive unless at least one of £3 or bi is nonzero. This is 
also the case for any 62 including 62 7^ 0, as remarked in Section 3] Nonzero £3 
or bi however leads to a weak singularity at z = 0. 



17 



We comment that our ' (z) is a generalization of the result of Celerier [9] 
which corresponds to the case = 0. To confirm the consistency, however, there 
are two remarks to make; (i) Ref. [9] employs the FLRW gauge, which is different 
from ours. To compare, an appropriate transformation is needed between the 
differential coefficients of the three LTB functions, (ii) Eq.(45) of [9] contains 
two wrong numerical factors; 7 and 10 should read, respectively, 6 and 9. We 
have confirmed that our result is consistent with [9] after taking these two points 
into account. 
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